Champion Statistics

This notebook outlines my process of generating models regarding champion statistics. This notebook is dependent on the data table gameInfo generated from DataExtraction.RMD.

Packages

library(plotly)

Warning: package ‘plotly’ was built under R version 4.1.2
Loading required package: ggplot2

Attaching package: ‘plotly’

The following object is masked from ‘package:ggplot2’:

    last_plot

The following object is masked from ‘package:stats’:

    filter

The following object is masked from ‘package:graphics’:

    layout

library(RColorBrewer)

Warning: package ‘RColorBrewer’ was built under R version 4.1.1

library(data.table)

data.table 1.14.0 using 4 threads (see ?getDTthreads).  Latest news: r-datatable.com

library(cluster)

Changing Data Structure and Making a List to Store Results

gameInfo.win <- gameInfo %>% 
  mutate(win = as.logical(win)) # Changing to logical to make future code easier
# # Commenting out to make sure that I don't run over my cashed results
# winrate <- list(
#   tables = list(),
#   plots = list()
# )

Overall Champion Winrate

winrate$tables$base <- gameInfo.win %>% 
  group_by(championName) %>% 
  summarize(winrate = mean(win), games = n(), .groups = "drop")

head(winrate$tables$base)

Lux is very popular and has a high win rate. Also yeesh, Ryze has a really low win rate in Season 11. Maybe we can see a bit more when we include the tier of play. Weirdly a lot of marksmen in the very popular category. From inspection there might be three groups of champions: 1. Balanced with medium play rate 2. Balanced with high play rate 3. Under tuned with low play rate My worry with a k-means is that the relative importance games played will overshadow the winrate statistic - maybe we can normalize the dataset.

Adding Tier Information

winrate$tables$tier <- gameInfo.win %>% 
  group_by(championName, tier) %>% 
  summarize(winrate = mean(win), games = n(), .groups = "drop")

head(winrate$tables$tier)

winrate$plots$tier <- winrate$table$tier %>% 
  plot_ly(
    x = ~games,
    y = ~winrate,
    color = ~tier,
    text = ~championName,
    type = "scatter",
    mode = "markers"
  )

winrate$plots$tier

Well we can for sure see regression towards the mean with a higher number of games wit ha few notable exceptions: 1. Rell is popping off in Diamond but likely variance due to low game number. Similar with Ornn in Iron. 2. One notable outlier that can be viewed instantly is Iron Yuumi, with a winrate of 0.434 which matches intuition. 3. Similarly, Ryze, Zoe, and Gwen have horrible win rates in bronze for the number of games. However this also might be due to oversampling of particular players in the scraping process.

A few ADC’s also grab my attention - Ezreal and Lucian seem to have an especially low winrate in Diamond. Might be worth investigating. This data set doesn’t have any obvious clusters to my eyes, let’s try running a kmeans clustering algorithm on the overall play rates.

K-means Clustering

Normalize Data Set and run kmeans models

It would appear that 3 or 4 clusters is optimal - let’s plot them: ## Plotting ### Setting up Data Table

Graphing Method 1

# Will use this in the future, even if there is that ugly play button
winrate$tables$clustered %>% 
  pivot_longer(cols = 4:11, names_to = "n_clusters", values_to = "membership") %>%
  plot_ly(
    x = ~games,
    y = ~winrate,
    colors = "Set3",
    color = ~membership,
    frame = ~n_clusters,
    type = "scatter",
    mode = "markers",
    text = ~championName
  ) %>% 
  layout(
    title = "Cluster Membership",
    xaxis = list(title = "Games Played"),
    yaxis = list(title = "Winrate")
  ) %>% 
  animation_opts(
    frame = 100
  )

Graphing Method 2 (Really annoying but better plot)

The 3 cluster model was too general and the 4 cluster model identifies an under performing group but includes champions like Cassiopiea with a decent win rate (0.495) in this group. The 5 cluster model appears to capture this under performing group of champions without gross overestimating. This is presuming that there is indeed an underlying structure to this space which may not be the case. It may be the case that a more complex space may yield more representative results as k-means clustering uses euclidean distance.

It’s also unclear whether or not the low play rate is a cause of the low win rate or because of champion imbalance. Optimal prescriptive balance changes might simply be a simple binary classifier in which champions below a certain win rate need buffs and the converse for high win rate champions.

So low play rate might be caused by a few factors: 1. The champion is weak and thus people avoid playing them. Champion weakness can be because: A. The champion cannot carry even when given gold B. They have a hard time acquiring gold Note that these factors could be cause by external factors like the metacgame or item changes. 2. The champion is difficult or un-fun and thus fewer players choose to play said champion.

To test hypothesis 1, we can control for gold and examine the winrates of the lowest winrate champions. If the champion’s winrate is no different from the null hypothesis (50%) at higher gold values, then sub hypothesis A is false. From there, we can look at the average gold distribution for that champion, and if it is significantly lower than the null hypothesis, then hypothesis B is true.

If both hypotheses are false, then we can use either internal or external (Mobalytics) measures of difficulty to determine the accuracy of hypothesis 2.

Preparation

High / Low WR table of champions

Let’s take the 3 champions per role with the lowest winrates and compare them to the top 3 winrate champions by role.

Function To Execute Model Generation

GET_WR_DATA <- function(ROLE, TABLE1 = data.temp$goldPercent, TABLE2 = winrate$tables$highlow){
  
  glm_data <- TABLE1 %>% 
    filter(
      championName %in% (TABLE2 %>% 
        filter(individualPosition == ROLE) %>% 
        {.$championName})
    ) %>% 
    return()
  
}

GET_WR_MODEL <- function(table){
  
  glm(win ~ championName*goldEarned, data = table, family = "binomial") %>% 
    return()
  
}

GET_WR_PLOTS <- function(ROLE, table, table2 = winrate$tables$highlow, ...){
  
  models <- list()
  preds <- list()
  goldRange <- data.frame(goldEarned = seq(0, .75, length.out = 75))
  
  champions <- table %>% 
    distinct(championName) %>% 
    left_join(
      table2, 
      by = "championName"
    )

  
  plot <- plot_ly()
  
  for(i in 1:nrow(champions)){
    
    models[[i]] <- table %>% 
      filter(championName == champions$championName[[i]]) %>% 
      filter(individualPosition == ROLE) %>% 
      glm(formula = win~goldEarned, family = "binomial", data = .)
    
    preds[[i]] <- predict.glm(models[[i]], type = "response", newdata = goldRange)
    
    if(champions$win[[i]] == "high"){
      
    plot <- plot %>% 
      add_trace(
        x = goldRange$goldEarned,
        y = preds[[i]],
        type = "scatter",
        mode = "lines",
        name = champions$championName[[i]]
      )
    } else if(champions$win[[i]] == "low"){
      
      plot <- plot %>% 
      add_trace(
        x = goldRange$goldEarned,
        y = preds[[i]],
        type = "scatter",
        mode = "lines",
        name = champions$championName[[i]],
        line = list(dash = "dash")
      )
      
    }
    
  }
  
  plot <- plot %>% 
    layout(
      xaxis = list(title = "Gold Percent"),
      yaxis = list(title ="Predicted Probability of Winning"),
      title = paste0("Top and Bottom 3 Winrate Champions (", str_to_title(ROLE), ")")
    )
  
  return(plot)
}

Observing Winrate when Controlling for Gold Percent

winrate$tables$WR_highlow$plot

[[1]]

[[2]]

[[3]]

[[4]]

[[5]]

NA

Covariate Champion Winrate

Are there any particular combinations of champions which are truly degenerate (see Master Yi - Taric funneling which was a gigantic issue in previous seasons) ## Temporary Table to Make Funciton More Efficient

data.temp$champTeams <- gameInfo.win %>% 
  select(match, win, championName) %>%
  group_by(match, win) %>% 
  mutate(champion = row_number()) %>% 
  pivot_wider(
    names_from = champion,
    values_from = championName
  ) %>% 
  mutate(team = str_c("_",`2`,`3`,`4`,`5`,"_", sep = "_")) %>% 
  select(match, win, champion = `1`, team) %>% 
  ungroup()

Executing Function

winrate$tables$covariate.temp <- expand.grid(
  champions.scraped$name,
  champions.scraped$name
) %>%
  apply(., 1, sort) %>% # Removing rows symmetrically, no need to double search
  t() %>% 
  unique() %>% 
  as_tibble() %>% 
  rename(Champ1 = V1, Champ2 = V2) %>% 
  filter(Champ1 != Champ2) # Removing diagonal elements

Warning: The `x` argument of `as_tibble.matrix()` must have unique column names if `.name_repair` is omitted as of tibble 2.0.0.
Using compatibility `.name_repair`.
This warning is displayed once every 8 hours.
Call `lifecycle::last_warnings()` to see where this warning was generated.

# This takes SO damn long to run, I wonder if there is a better way to do this...
# Doing this in two steps to avoid running long functions again
winrate$tables$covariate.temp2 <- winrate$tables$covariate.temp %>%
  mutate(
    gameList = pmap(
      .,
      .f = function(Champ1, Champ2, TABLE = data.temp$champTeams){

        TABLE %>%
          group_by(match, win) %>%
          filter(champion == Champ1) %>% # Okay this is way faster but still turbo slow
          filter(str_detect(team, paste0("_", Champ2, "_"))) %>%
          ungroup() %>%
          return()

      }
    )
  )

winrate$tables$covariate <- winrate$tables$covariate.temp2 %>% 
  mutate(
    data = map(
      gameList,
      .f = function(gameList){
        
        gameList %>% 
          summarize(winrate = mean(win), games = n()) %>% 
          return()
        
      }
    )
  ) %>% 
  unnest(cols = c(data))

# Making a new table to make searching particular champion pairs easier
winrate$tables$covariate.search <- bind_rows(
  winrate$tables$covariate,
  winrate$tables$covariate %>% 
    rename(Champ2 = 1, Champ1 = 2)
)


winrate$plots$covariate <- winrate$tables$covariate %>% 
  plot_ly(
    x = ~games,
    y = ~winrate,
    text = ~paste(Champ1, Champ2, sep = " "),
    type = "scatter",
    mode = "markers"
  )

winrate$plots$covariate 

Warning: Ignoring 2510 observations
Warning: Ignoring 2510 observations

Probably not enough games to make any predictions - with 155*154 champion duo combinations, would need a lot more data to have a decent sample of all combinations.